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Effect of tube diameter and capillary number on platelet 
margination and near-wall dynamics 

Timm Kruger 


Abstract The effect of tube diameter D and capil¬ 
lary number Ca on platelet margination in blood flow 
at ~ 37% tube haematocrit is investigated. The sys¬ 
tem is modelled as three-dimensional suspension of de¬ 
formable red blood cells and nearly rigid platelets us¬ 
ing a combination of the lattice-Boltzmann, immersed 
boundary and finite element methods. Results show 
that margination is facilitated by a non-diffusive ra¬ 
dial platelet transport. This effect is important near 
the edge of the cell-free layer, but it is only observed for 
Ca > 0.2, when red blood cells are tank-treading rather 
than tumbling. It is also shown that platelet trapping 
in the cell-free layer is reversible for Ca < 0.2. Only for 
the smallest investigated tube {D = 10 pm) margination 
is essentially independent of Ca. Once platelets have 
reached the cell-free layer, they tend to slide rather than 
tumble. The tumbling rate is essentially independent of 
Ca but increases with D. Tumbling is suppressed by the 
strong confinement due to the relatively small cell-free 
layer thickness at ~ 37% tube haematocrit. 

Keywords Platelet margination • red blood cell • 
cell-free layer • lattice-Boltzmann method • immersed 
boundary method • simulation 
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1 Introduction 

The aim of this article is to investigate the effect of 
tube diameter and capillary number (wall shear rate) on 
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the margination and near-wall dynamics of platelets in 
blood flow via highly resolved three-dimensional com¬ 
puter simulations. Platelet adhesion or activation are 
not in the focus of this research. 

Human blood is, by volume, composed of about 55% 
plasma (mostly water and proteins) and 45% suspended 
blood cells. The majority of the blood cells are red blood 
cells (RBCs), also called erythrocytes. Other, less com¬ 
mon blood cells are white blood cells (leukocytes, one 
for 1000 RBCs) which form a significant part of the 
system ( [Robertson et al[ 2007|) and platelets 


immune 


(thrombocytes, one for 15 RBCs) (AlMomani et al 


2008| ) . The volume fraction of the cellular phase is called 
haematocrit Ht. 

Platelets play an important role in the process of 
blood clotting and the repair of damaged vessel walls 
(endothelium). By mechanical contact with damaged 
endothelium, platelets can be activated. The activation 
involves a reorganisation of the platelet cytoskeleton. 


making them more flexible (Fogelson and Guy 2008). 


Activated platelets are capable of forming networks via 
fibrinogen and von Willebrand factor molecules giv¬ 


ing rise to the build-up of blood clots (Kulkarni et al 


2000 Doggett et al 2002). Non-activated platelets can 


be considered as rigid discoid particles over the en¬ 
tire physiological shear stress range (up to a few 10 Pa) 


(Goldsmith and Mason 1967 Teirlinck et al 1984[ |Tn-] 


ritto and Goldsmith 1992). 


The main task of platelets is to recognise regions 
of damaged endothelium and to assist in its repair. In¬ 
deed, it is observed that the platelet concentration is, 
under certain conditions, significantly increased close to 


the walls (Turitto et al 1972 Goldsmith and Turitto 


1986). This near-wall excess is of physiological impor¬ 


tance as it increases the probability for the platelets 


to adhere to the damaged endothelium (Yeh and Eck¬ 


stein 


1994). Platelet concentration near the wall can 


be several times higher than that at the centre of the 
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flow ( Tangelder et al[ 1982[ Tilles and Eckstein 1987 


Aarts et al[ 1988 Eckstein et al[ 1989). The lateral 
motion of platelets, leading to this concentration in¬ 
homogeneity, is called margination. Usually, the near¬ 
wall concentration peak is located within or near the 


cell-free layer (CEL) which is depleted of RBCs (Yeh 
and Eckstein! p!9M| ) . The CEL thickness ^cfl plays an 
important role in blood flow, leading to the Eahraeus 


and Eahraeus-Lindqvist effects (Eahraeus and Lindqvist 


1931 Lei et al 2013). 


RBCs have a twofold effect on the platelet motion. 
Eirst, it has been observed that the presence of RBCs 
increases the diffusivity of platelets in blood vessels by 
orders of magnitude compared to their Brownian diffu¬ 


sivity (Aarts et al 1983 1984 Breugel et al 1992 Tu- 
ritto and Goldsmith| 1992). Eurthermore, non-diffusive 
platelet margination towards the vessel wall increases 


strongly when RBCs are present (Turitto and Weiss 
198Q| Cadroy and Hanson] 199Q| [Joist et al[ 1998||Peer- 


schke et al 


2004) 


The wall shear stress and the haematocrit have been 


identified as important parameters for margination (Tu- 

ritto and Baumgartner 

1975[ |1979| Turitto and Gold- 

smith| |1992| [Zhao et al 

|2007| [Eogelson and Guy| 2008 

Reasor Jr et al 2013). 

It is known that margination is 


only significant when the wall shear rate is above 200 s 


-1 


(Tilles and Eckstein 1987 Eckstein et al 1988 Bilsker 


et al[ 1989). There also seems to be a minimum haemat 
ocrit of 10% below which no margination is observed for 


physiological stresses (Tilles and Eckstein 1987 Waters 


and Eckstein 1990). Margination becomes more impor¬ 


tant with increasing haematocrit (Yeh and Eckstein 


1994). Also platelet size and shape play a role for the 


margination efficiency (Eckstein et al 1988 Thompson 
et al 2013 Reasor Jr et al[ 2013 Muller et al[ 2014). 

Only in recent years, resolved computer simula¬ 
tions became available to study the mechanisms of 
platelet margination directly, without assuming effec¬ 
tive platelet transport models. Simulations are able to 
provide microscopic details (such as stresses, cell ve¬ 
locities and deformations) that are difficult to measure 
in experiments. As already assumed in the 1970s by 


the RBCs are sufficiently soft, these collisions lead to 
directed events which are necessary to explain margina¬ 


tion (Eckstein and Belgacem 1991). It can be concluded 


that the dynamical state of RBCs (tumbling versus 
tank-treading) plays a major role in platelet margina- 


tion, as already suggested by Yeh and Eckstein (1994). 

Recent review articles provide more details ( 

Kumar and 

Graham 2012a Eogelson and Neeves 

2015 

). 


In general, platelet margination is still not well un¬ 
derstood. This research aims at investigating in more 
detail the effect of tube diameter and capillary num¬ 
ber on margination. Three-dimensional computer sim¬ 
ulations employing the lattice-Boltzmann, immersed 
boundary and finite element methods (section are 
performed to investigate platelet trajectories in an en¬ 
vironment of deformable RBCs. The relevant parameter 
space and the system setup are presented in section 
The results about platelet margination and platelet dy¬ 
namics in the CEL are presented and discussed in sec¬ 
tion 1^ Einally, the observations and implications are 
summarised in section [H 


2 Physical model and numerical methods 

The numerical methods employed in this work are 
briefly outlined. Details and benchmark tests are re- 

Kriiger et 


ported elsewhere (Kruger 2011 


2013, 2014). The lattice-Boltzmann method (section 


2.1) is used as Navier-Stokes solver, and the immersed 
boundary method (section [2^ assumes the role of the 
fluid-structure interaction. Eor the membranes of the 
red blood cells, a finite element approach (section [2.3| ) 
is employed. In section [2^ the numerical implementa¬ 
tion of the viscosity contrast is outlined. Tab. [^provides 
an overview of all relevant symbols and parameters. 


2.1 Lattice-Boltzmann method 

The standard lattice-Boltzmann method with D3Q19 
lattice ( Qian et~aH 1992) and BGK collision operator 
( Bhatnagar et al[ 1954) is employed to solve the Navier- 


Goldsmithl (fTMl) and iTuritto et alj ([19^, it is nowa- Stoke s equations ( |Succi[ |2QQ1| [Aidun and Clausen 


days generally accepted that the dynamics of RBCs and 
their hydrodynamic interactions with platelets under 
the influence of shear are responsible for margination. 


2010 ). 


rather than a volume exclusion effect (Crowl and Eogel 


son 


2010). Zhao and Shaqfeh (2011) found that the lat 


The BGK operator is parametrised by a relaxation 
time xAt (with At being the time step). The kinematic 
fluid viscosity is directly related to T according to 


eral platelet motion in the RBG-rich region is diffusive 
and caused by shear-induced fluid velocity fluctuations. 
According to Kumar and Graham (2012b), platelet 


1 / \\Ax^ 

^“3 V^~2 ) 'aT 


( 1 ) 


margination, at least in the dilute regime, is caused 
by heterogeneous collisions of platelets and RBGs. If 


where Ax is the lattice resolution. 

External forces, such as those from the immersed 


boundary method (section 2.2), are coupled through 
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(a) (b) (c) 


Fig. 1 Mesh generation for a red blood cell. The 20 face ele¬ 
ments of an icosahedron (a) are subdivided, and the resulting 
new vertices shifted to the circumsphere (b). The mesh of a 
red blood cell (c) can be obtained by adjusting the vertex po¬ 
sitions in space. The meshes in (b) and (c) have 2880 elements 
each. 


the forcing scheme by Shan and Chen (Shan and Chen 


1993). 


2.2 Immersed boundary method 


The immersed boundary method (IBM) ( [Peskin 1972 


2002) is used to couple the fluid flow on the Eulerian 


lattice and the off-lattice membrane interface dynam¬ 
ics. IBM has often been employed for the simulation of 


elastic deformable objects, e.g. by Eggleton and Popel 


(1998); Zhang et al (2007 ); |Sui et al (2008); Eogelson 
and Guy (2008); Doddi and Bagchi (2009). 


Eluid velocities are interpolated at the locations 
of the membrane vertices, and membrane forces (sec¬ 


tion 2.3) are distributed back to the lattice. Eor this 


purpose, a short-range tri-linear interpolation stencil is 
used that requires only 2^ lattice points in three dimen¬ 


sions (Peskin 2002) 


2.3 Red blood cell mesh and elasticity model 


In their undeformed state, RBCs are biconcave discs 
with a diameter of about 8 pm. RBC membranes are 
characterised by two major elastic contributions: in¬ 


plane shear and bending resistance (Skalak et al 


1973 


Evans 


1974[ ). Eurthermore, local and global area 


changes are small (typically below 1%) since the lipid 
bilayer is incompressible. The total RBC surface ^ 
135 pm^ is therefore essentially constant. Under phys¬ 
iological conditions, the RBC volume, ?^94pm^. 


is also constant (Evans and Eung 1972). The strong 


deformability of RBCs is facilitated by the large area 
excess (additional area compared to a sphere with the 
same volume) which is about 35%. This way, RBCs 
can squeeze through capillaries with diameters as small 
as 4 pm without violating the surface and volume con¬ 


straints (Skalak and Branemark 1969) 


In the present work, RBCs and platelets are treated 
as closed membranes discretised by Nf flat faces (area 
elements). Each element consists of three nodes (ver¬ 
tices) of which there are in total A^n = (A^f + 4)/2. The 
generation of the RBC mesh, depicted in Eig. is 
based on the icosahedron-refinement procedure (Ra¬ 


manujan and Pozrikidis, 1998 Kruger 2011). In the 


present work, each RBC mesh consists of 2880 elements 
and 1442 vertices. Platelets are assumed to be non- 
activated and modelled as nearly rigid, oblate ellipsoids 
with an aspect ratio 0.28 and major radius of 1.8 pm. 
Each platelet mesh comprises 320 faces and 162 ver¬ 
tices. 

The RBC membrane is modelled as a hyperelastic 
continuum with four elastic energy contributions as out¬ 
lined below. Membrane viscosity is neglected. It is as¬ 
sumed that undeformed RBCs are stress-free. The RBC 
equilibrium shape is imposed through the input geom¬ 
etry in Eig. |l(c)| 

Skalak’s constitutive law for the in-plane energy 


density is employed in the present simulations (Skalak 


et al 1973): 


Es = I{if+ 2h- 2h) + ^/|) . 


( 2 ) 


The two parameters and Ka control the strength 
of the membrane response to local shear deformation 
and dilation, respectively. Eor healthy R BCs, the values 
are Ks =5.3pNm“^ and = 0.5Nm“^ (Gompper and 


Schick 2008). The in-plane deformation is described 


by the two strain invariants: Ii and I 2 . Eq. ([^ is discre¬ 
tised on the membrane mesh according to |Charrier et al 
(1989) and Shrivastava and Tang (1993). More details 
are provided by Kruger et al (2011). 

Eor the bending energy, the discretised form 


Eb = 


V3kb 




( 3 ) 


{ij) 


is used. Here, Kb is the bending resistance (2-10 N m 
for an RBC), Gij is the angle between the normal vec¬ 
tors 11 / and hj of two neighbouring elements, and is 
the corresponding equilibrium angle of the undeformed 
membrane. The sum runs over all pairs of neighbouring 
elements. 

In order to maintain a nearly constant RBC surface 


area and following Evans and Skalak (1980) and Seifert 


(1997), an additional energy is introduced: 


Ea = 


Ka 


(a-A(o))' 


A(0) 


( 4 ) 


The magnitude of the surface energy is controlled by 
the surface modulus k*a. A similar approach is used for 
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Fig. 2 Indicator field /(x) for distinguishing between inte¬ 
rior and exterior fluid sites. A two-dimensional cross-section 
of a three-dimensional RBC suspension in a tube is shown. 
Red lines denote the Lagrangian membranes. The grey-scale 
lattice-like portion of the figure reflects /(x) on the same cross- 
section. Black denotes the interior, white the exterior. 


is the viscosity contrast. Fig.j^shows the indicator field 
for a given RBC configuration. 


3 System parameters and simulation setup 


The system parameters and the simulation setup are 
discussed. After introducing the relevant system param¬ 
eters in section [3T| the choice of simulation parameters 


is briefly discussed in section |3.2| The initialisation of 
the simulations (section 3.3) deserves some additional 
attention. 


3.1 Physical parameters 


the RBC volume (Seifert 1997): 



(5) 


The total membrane energy E is the sum of all con¬ 
tributions Fs, Fb, Fa and Fy. The forces acting on 
membrane vertices at position Yj are computed from 
the energy functional F({rj}) via the principle of vir¬ 
tual work: 


, a£({r,}) 

^ drj 


(6) 


The differentiation is performed analytically, and the 
resulting forces are implemented in the code. A detailed 
derivation of all force contributions is given by |Krhger| 

dMT]). 


2.4 Viscosity contrast 


In order to implement a viscosity contrast A, i.e. dif¬ 
ferent viscosities inside and outside of the RBCs, each 
lattice site requires up-to-date information about its lo¬ 
cation relative to nearby Lagrangian vertices. A signed 


distance field as detailed by Frijters et al (2015) is used 


to construct an indicator field /(x) at each lattice site x, 
denoting the relative location of each site with respect 
to nearby RBCs and platelets (7 = 0 outside and 7=1 
inside particles with a linear slope across the particle 
surface). The instantaneous kinematic viscosity is then 
obtained from the indicator field 


v(x) = Vout (1 -/(x)) + Vin/(x) 


(7) 


where Vout and Vin are the external and internal vis 
cosities, respectively, and 

Vin 


A = 


Vout 


( 8 ) 


The present research concerns with blood flow in 
straight tubes, which involves a series of parameters 
and related symbolds as collected in Tab. 

The flow is driven along the tube axis in positive 
v-direction by a pressure gradient p', mimicked by a 
constant and homogeneous force density / = In the 
absence of particles, this would lead to a parabolic ve¬ 
locity profile with peak velocity uq = p'7)^/(16T7ex), av¬ 
erage velocity uq = uqI2 and wall shear rate TV = ^uq/D. 
RBCs in such a flow experience viscous shear stresses 
that can be characterised by the capillary number 

_ Ve^fw^RBC _ p'D^rbc /gx 

Ks ~ 4 ks ' 


Table 1 Overview of used parameters and symbols. 


Parameter 

symbol 

value 

Tube radius, diameter 

R, D 

variable 

Tube segment length 

L 

48 pm 

RBC radius 

rRBC 

4 \im 

Platelet radius 


1.8 pm 

Platelet thickness 

^pl 

1 pm 

Tube haematocrit 

Htt 

^37% 

RBC count 

N-rbC 

variable 

Platelet count 

A^pl 

Npi =Nbbc/2 

Fluid density 

P 

1000 kg m“^ 

Cytoplasm viscosity 

flin 

5mPas 

Plasma viscosity 

flex 

1 mPas 

viscosity contrast 

A 

5 

RBC shear elasticity 

K-s 

5.3pNm-i 

RBC bending rigidity 

Kb 

2-10-i9Nm 

Reduced bending modulus 

KB/{Ksr^) 

1/424 

Pressure gradient 

p' 

variable 

Capillary number 

Ca 

variable 

Lattice resolution 

Ax 

0.33 pm 

Time step 

At 

variable 

Centre velocity (no cells) 

Uq 

().05Ax/At 

Average velocity (no cells) 

Uq 

uqI2 

Wall shear rate 

Yw 

Auq/D 

Number of time steps 

Nt 

2-105 

Advection time 

^ad 

96()At 
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Under physiological conditions (T]ex, K's and tkbc as re¬ 
ported in Tab.[^, Ca=l corresponds to 7w ~ 1300s“\ 
which leads to significant RBC deformations. Note that 
the definition of Ca in Poiseuille flow is not unique and 
that there may be a constant conversion factor between 
the present and other authors’ choice. 

The purpose of this article is to investigate the effect 
of tube diameter D and capillary number Ca on the 
platelet dynamics in tube flow at physiological blood 
parameters. Tab. lists the investigated values for D 
and Ca. The tube haematocrit is kept at Htt ^37% in 
all cases. 

In order to And a dimensionless time for the advec- 
tion of the suspension, one may define the advection 
time scale 


^ad — 


2rRBC 

uo 


( 10 ) 


It is the average time an RBC requires to move its own 
diameter in the unperturbed flow. 


3.2 Simulation parameters 

The flow is periodic along the tube axis (v-direction). 
The circular tube cross-section is approximated by a 


staircase, and the half-way bounce-back boundary con¬ 
dition is used to enfore no-slip at the wall (Ladd 1994). 
Due to the large tube diameters and the choice of the 
external viscosity, artificial numerical slip is negligible. 

The lattice resolution is Ax = 0.33 pm, i.e. the RBC 
diameter is 2rBBC = 8 pm = 24Ax. While the tube diam¬ 
eter D varies, the length of the simulated tube segment 
is L = 12rBBC = ^44Ax in all simulations. 

Obeying A = 5, the viscosities are Vex = 1/6 and 
Vin = 5/6 in simulation units. This leads to BGK re¬ 
laxation parameters Tex = 1 and Tin = 3, according to 


Eq. 

The flow is driven by a constant force density f = p' 
along the tube axis. Its value is chosen in such a way 
that the centre velocity uq would be 0.05 Av/A^ in the 
absence of any particles. This is to avoid compressibility 
effects and to keep the time step sufficiently small to 
achieve stable simulations. Note that the actual centre 
velocity is smaller (typically by a factor of ^ 2) due to 
the presence of the cells. 

Tab. lists all relevant parameter values for the sim¬ 
ulations undertaken. 


3.3 Simulation initialisation 


Table 2 Overview of numerical simulation parameters, [s.u.] 
denotes “simulation units”. Symbols are explained in Tab. 
Other simulation parameters that are the same throughout 
are Vin =5/6, Vext = 1/6? ~ 37%, = 0.5, /Ca = Xy = 1 and 

A/, =2-10L 


D 

Ca 

Npi 

p' 

KS 


[Hui] 



[s.u.] 

[s.u.] 

[s.u.] 


0.1 



0.133 

0.0453 


0.2 



0.0666 

0.0226 

10 

0.3 

7 

1.48-10-4 

0.0444 

0.0151 

0.6 

/ 

0.0222 

0.00755 


1.0 



0.00453 

0.00453 


2.0 



0.00666 

0.00226 


0.1 



0.0888 

0.0302 


0.2 



0.0444 

0.0151 

15 

0.3 

16 

6.58-10-5 

0.0296 

0.0101 

0.6 

0.0148 

0.00503 


1.0 



0.00888 

0.00302 


2.0 



0.00444 

0.00151 


0.1 



0.0666 

0.0226 


0.2 



0.0333 

0.0113 

20 

0.3 

28 

3.70-10-5 

0.0222 

0.00755 

0.6 

0.0111 

0.00377 


1.0 



0.00666 

0.00226 


2.0 



0.00333 

0.00113 


0.1 



0.0444 

0.0151 


0.2 



0.0222 

0.00755 

30 

0.3 

63 

1.65-10-5 

0.0148 

0.00503 

0.6 

0.00741 

0.00252 


1.0 



0.00444 

0.00151 


2.0 



0.00222 

0.000755 


Due to the relatively small system size (between 
Nrbc = 14 and 126), the number of platelets has been 
taken as N^\ =Nrbc/2. This is about 7-8 times larger 
than observed under physiological conditions. Since 
platelets are much smaller than RBCs, their volume 
fraction is still small compared to Htt. 

At the beginning of a simulation, platelets and 
RBCs are distributed throughout the tube with ran¬ 
dom positions and orientations. Any overlap of parti¬ 
cles with particles and particles with the tube wall is 
avoided. The platelets are positioned first, which as¬ 
sures a more homogeneous distribution of them across 
the tube cross-section. In order to facilitate this proce¬ 
dure, all particles are initially shrunk to half of their lin¬ 
ear size and grown afterwards within 4000 time steps. 
During growth, the particle volume is increased with 
a constant rate and overlap is avoided by a repulsion 
force. 

As the growth process leads to an increase of kinetic 
energy in the system, a friction force is used to dissi¬ 
pate energy. This force is switched off once the growth 
process is complete. The exact form of the dissipation 
force is not significant for this purpose. 

After all particles have reached their full size, the 
simulation starts with a constant fluid density p (unity 
in simulation units) and zero velocity everywhere. The 
force density f = p' along the tube is switched on in- 
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stantaneously to drive the flow. Each simulation runs 
for Nt =2' 10^ time steps, i.e. about 208rad- 


4 Simulation results and discussion 


The results and discussion are presented in three parts. 


In section [4T| the final suspension state is visually in¬ 
spected and characterised, and the cell-free layer thick¬ 
ness is investigated. Section [4^ contains the margina- 
tion observations and analysis. Platelet dynamics near 
the tube wall is discussed in section [T3l 


4.1 Visual results and cell-free layer thickness 

Fig. i shows the final state (after 208 advection times) 
for the smallest and largest simulated capillary num¬ 
ber Ca and each tube diameter D. Both the RBC and 
platelet configurations depend on the geometry and 
flow parameters. 

At Cn = 0.1, RBCs are only weakly deformed, and 
their dimples are easy to spot. Contrarily, at Ca = 2.0, 
RBCs show a strong and characteristic elongation, and 
the dimples have disappeared. The CEL is particulary 
pronounced at higher Cn, and RBCs close to the tube 
wall tend to align with the wall and form a ring. This 
leads to a slight increase in the local RBC concentra¬ 
tion (Fig. |4(a)] ) as observed previously (Lei et al 2013). 
Furthermore, RBCs in Fig. |3(b)| show a typical zigzag 
configuration that can be found in small blood vessels 


(McWhirter et al 2009). 


The platelet distribution changes strongly with D 
and Cn, as discussed in more detail in section |4.2[ In 
summary, platelets prefer to marginate towards the 
tube wall with increasing Ca. Also, for smaller D, a 
larger fraction of platelets resides near the tube wall at 
the end of the simulations. 

There is no unique way to define the CFL thick¬ 
ness ^CFL since the edge between the outermost RBCs 


and the plasma layer is diffuse (Lei et al 2013). The 
blurring of the CFL is clearly visible in Fig. In the 
following, the definition of ^cfl is based on the radial 
tube haematocrit profile averaged over time. As 

shown in Fig. 4(a)[ Htt{r) decreases strongly near the 
wall. The CFL thickness is defined by the distance from 
the tube wall where reaches Htt/2. 

Fig. |4(b)| shows the resulting CFL thickness as func¬ 
tion of Ca for different tube diameters D. £cfl is 
generally increasing with Ca until it saturates above 
Ca ~ 0.3-0.6. This is close to the point where RBC 


tumbling is replaced by tank-treading (Kruger et al 


to their large deformations. Evidently, the CFL thick¬ 
ness is larger for wider tubes, but the relative thick¬ 
ness icFh/D is decreasing. In the macroscopic limit, for 
D > 100 pm, the CFL is less important for the effective 


blood rheology (Lei et al 2013). 


2013). The slight decrease of ^cfl at Cn = 2 is caused 


In the following, ^CFL is used to identify those 
platelets that are located between the RBC-rich region 
and the tube wall. 


4.2 Platelet margination 

Some exemplary radial platelet positions as function of 
time are shown in Fig. In the following, a platelet is 
considered being fully marginated when it is located in 
the CFL vicinity with its radial centre position between 
/? —2^CFL and R (indicated by dotted lines in Fig. |^. 
As seen in Fig. |4(a)[ this corresponds to the region be¬ 
tween the tube wall and the haematocrit peak at the 
CFL edge. In this region, Htt{r) is a strongly decreasing 
function. 

The character of the platelet motion depends on the 
tube diameter D and the capillary number Cn. |Zhao and| 
Shaqfeh (2011) demonstrated previously that platelets 
show diffusive behaviour in the RBC-rich region until 
they reach the CFL where they get trapped irreversibly. 
However, the present data shows that this is only the 
case for larger capillary numbers {Ca > 0.2). Fig. [^re¬ 
veals that, for Cn = 0.1, several platelets leave the CFL 
vicinity again. Those platelets do not move far away 
from the CFL, though; instead they tend to become 
caught again shortly after. For larger Cn, platelets are 
trapped in the CFL once they visit it for the first time. 

Since it is difficult to extract further information 
directly from the trajectories in Fig.|^ Fig.j^shows the 
time evolution of the average radial platelet position 
and the fraction of platelets that are located within 
2 ^cfl from the tube wall. The first observation is that 
the time for platelets to become trapped increases with 
tube diameter. In larger tubes, platelets have to move 
for a longer distance to reach the CFL; this takes more 
time. From the third and fourth rows in Fig. one can 
conclude that margination is still ongoing for D>20 pm 
at the end of the simulation (after 208 advection times). 
On the contrary, the platelet distributions in the smaller 
tubes with D < 15 pm seem to have reached a quasi¬ 
equilibrium. 

Secondly, margination and trapping are faster for 
higher capillary number. However, there is no strong 
difference between Ca = 0.3 and Ca = 2.0 as the second 
and third columns in Fig. show. Zhao and Shaqf^ 
(2011) observed that the platelet diffusion in the RBC- 


by a subtle geometrical rearrangement of the RBCs due rich region does not change significantly between Ca = 






































Effect of tube diameter and capillary number on platelet margination and near-wall dynamics 


7 



(a) D = 10 pm, Ca = 0.1 



(c) D =\5 pm, Ca = 0.1 



(e) D = 20pm, Ca = 0.1 



(g) 2) = 30 pm, Ca = 0.1 



(b) D = 10 pm, Ca = 2.0 



(d) D = 15 pm, Ca = 2.0 



(f) D = 20pm, Ca = 2.0 



(h) D = 30 pm, Ca = 2.0 


Fig. 3 Configuration of the red blood cells (red) and platelets (yellow) for all investigated tube diameters {D = 10, 15, 20, 
30 pm) and the smallest and largest studied capillary numbers (Ca = 0.1, 2.0) at Htt ~ 37% after t = 208 ^ad- The left columns 
show sections along the tube, the right columns tube cross-sections with reduced RBC opacity to reveal all platelets. Flow is 
from left to right (left columns) and out of the image plane (right columns). 


1.0 and Ca = 2.0, a fact that is attributed to the sat¬ 
uration of RBC deformation due to the conservation 
of the RBC surface area. The strong increase of the 
margination efficiency between Ca = O.l and Ca = 0.3 
in Fig. could be related to the dynamical state of 
the RBCs. RBCs tend to tumble below and tank-tread 
above Ca ~ 0.2 (Kruger et al 2013). Reasor Jr et al 


(2013) hypothesised that tank-treading RBCs may in¬ 
crease the margination efficiency. In fact, for Ca = 0.1 
(first column in Fig.[^, margination within the first 208 
advection times is nearly absent, except for the smallest 
tube with D = 10 pm. 

The smallest tube (D = 10 pm) plays a special role. 
Platelets tend to reach the CFL nearly instantaneously, 
even at Ca = 0.1. Therefore, margination in small ves¬ 
sels with D < 10 pm seems to be very effective for all 
capillary numbers. This is in line with the observa¬ 
tions by Reasor Jr et al (2013) who reported that, for 
rigid RBCs, margination is only observed for platelets 
that are not farther away from the wall than about 
5 pm. Therefore, all platelets in a lOpm-tube are al¬ 
ways in wall vicinity. This will be further discussed be¬ 
low. Already for a slightly larger tube (D = 15 pm), sev¬ 
eral platelets remain in the RBC-rich region for an ex¬ 


tended period of time if the RBCs are sufficiently rigid 
{Ca = 0.1). For tubes with D = 20 pm and 30 pm, only 
a fraction of the platelets initially located in the RBC- 
rich phase diffuse into the CFL during the simulation 
time. 

As a final remark regarding Fig. the number ra¬ 
tio of platelets near the tube wall fluctuates more for 
smaller Ca. This supports the observation that platelets 
are not irreversibly trapped in the CFL for Ca = 0.1. 
Further data substantiating this claim will be discussed 
shortly. 

Fig. shows some radial platelet distributions aver¬ 
aged over the final 25% of the simulation time (between 
~ 150 and 200 advection times). These data support 
the previous interpretations. For D = 10pm (first row), 
margination is very effective, and nearly all platelets are 
located close to the CFL (vertical dotted line). The situ¬ 
ation is different for larger tube diameters with Ca = 0.1 
(left column). In those cases, there is nearly no margina¬ 
tion, and the platelet distribution is close to being ho¬ 
mogeneous (horizontal dashed line). Much longer sim¬ 
ulations are necessary to decide whether margination 
for small Ca is just slow or whether this is already the 
equilibrium state. For Ca = 0.3 and Ca = 2.0 (second 
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(a) Radial haematocrit distribution {D = 30 p-m, Ca = 2.0) 



Ca 

(b) CFL thickness 


Fig. 4 Cell-free layer thickness ^CFL- (a) Exemplary radial tube haematocrit profile Htt{r) averaged over time. The CFL 
thickness ^CFL (denoted by grey area) is defined as the distance from the wall (r = R) to the point where Htt (r) reaches half of 
the average tube haematocrit. (b) ^CFL as function of capillary number Ca for all investigated tube diameters D at fftt ~ 37%. 
Error bars are estimated. Lines are guides for the eye. Note that the Ca-axis is logarithmic. 


and third columns in Fig. [^, platelets are effectively 
transported towards the tube wall. 

In Fig. [^h,i,k,l), a dip in the platelet concentra¬ 
tion close to the CFL edge is visible. This suggests 
that platelets are effectively transported in outward di¬ 
rection when they are close to the CFL edge, there¬ 
fore decreasing the concentration at r ^ R — 2^cfl and 
increasing it within the CFL. Platelet diffusion leads 
to a fresh supply of platelets moving from the RBC- 
rich region towards the wall until there are no platelets 
left in the tube interior. Further evidence supporting 
this interpretation will be discussed below. As men¬ 
tioned before, this process seems to have completed in 
Fig. me ,f), where the tube radius is relatively small 
(D = 15 pm). Yet, the process appears to be still ongo¬ 
ing in Fig. ,k,l), where the tube is larger. In the 
latter cases, it is expected that all platelets will even¬ 
tually reach the CFL after a sufficiently long time. 

As a remark, there is also a platelet concentration 
peak near the centreline in Fig.[^i,k). Note that the ab¬ 
solute number of platelets in that region is small since 
the bin area decreases linearly with r for r —)► 0. There¬ 
fore, the platelet concentration peak near the CFL edge 
is caused by many more platelets and is therefore more 
significant than the peak near r = 0. This is nicely borne 
out in Fig. |3(f)| where only three platelets are located 
near the centreline, but the remaining 25 are close to 
the CFL. Therefore, the central peak is probably a sta¬ 
tistical fluctuation due to the relatively short averaging 
period and limited platelet number. 


To investigate the radial platelet transport in more 
detail. Fig. shows the average radial platelet velocity 
(v^) and its standard deviation as function of radial 
position r for D = 20 pm and D = 30 pm. Due to the 
small number of platelets and the fast margination, the 
radial velocity data for D = 10 pm and D = 15 pm is too 
noisy to provide useful information. 


The first striking observation from Fig.j^is that the 
radial platelet velocity is close to zero everywhere, ex¬ 
cept near the edge of the CFL (grey area) for Ca > 0.2. 
This shows that platelets are pushed towards the wall if 
they have reached the edge of the CFL. This effect, how¬ 
ever, only exists when RBCs are sufficiently strongly 
deformed and in the tank-treading state. For the pre¬ 
sented tube diameters, collisions between platelets and 
nearly rigid RBCs {Ca < 0.2) do not lead to a significant 
lateral platelet transport near the CFL. Note that the 
lateral platelet velocity is closely related to the drift 


term in effective platelet transport models (Eckstein 
and Belgacem[ p^91| |Yeh and Eckstein 1994). 


The peak of (v^) (r) coincides with the peak of Htt{r) 
in Eig. |4(a)| Therefore, platelets are already pushed out¬ 
wards before they have passed the CEL edge. Eor both 
tube diameters shown in Eig.[^ (v^) (r) becomes signifi¬ 
cantly larger than zero when the distance to the wall is 
about 5 pm; this is in good agreement with the results 
report by |Reasor Jr et al (2013). 

Eig. [^b,d) reveals that the radial platelet velocity 
fluctuation strongly decreases within the CEL. This is 
an indicator for reduced platelet diffusion in that region. 
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Fig. 5 Platelet margination for different tube diameters {D = 10 pm with 7 platelets, D = 15 pm with 16 platelets, D = 20 pm 
with 28 platelets and P) = 30pm with 63 platelets). The radial position r (normalised by tube radius R) for each platelet is 
shown for two different capillary numbers {Ca = 0.\ and 2.0) as function of dimensionless time. The dashed line denotes the 
cell-free layer thickness with distance ^CFL from the wall. The dotted line denotes twice the CFL thickness, 2 ^cfl- Note that 
platelets are randomly distributed over the pipe cross-section at the start of the simulation. This means that more platelets 
are initially located at larger radii r since the cross-section area element obeys dA = 27rrdr. 
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t/t&d. 

(j) D = 30|J,m, Ca = 0.1 (k) Z) = 30|J.m, Ca = 0.3 (1) Z) = 30 |um, Ca = 2.0 

Fig. 6 Average platelet margination for all studied tube diameters D (varied row-wise) and a few capillary numbers (Ca = 0.1, 
0.3, 2.0, varied column-wise). Dashed lines show the radial position (r) averaged over all platelets as function of time. The grey 
regions bounded by solid lines indicate the number fraction {n) of all platelets that are located within a region 2 ^cfl from the 
tube wall (indicated by horizontal dotted lines). 
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r [|xm] 

(a) D = 10 pm, Ca = 0.1 



(d) D= 15 pm, Ca = 0.1 



(g) D = 20 pm, Ca = 0.1 



0 2 4 6 8 10 12 14 

r[pm] 

(j) D = 30pm, Ca = 0.1 



r [pm] 

(b) D = 10 pm, Ca = 0.3 



(e) D=15 pm, Ca = 0.3 



(h) D = 20 pm, Ca = 0.3 



r[pm] 

(k) 0 = 20 pm, Ca = 0.3 




1^ "1 

0 12 3 


4 5 

r[pm] 



(c) D =10 pm, Ca = 

2.0 




r[pm] 

(i) D = 20pm, Ca = 2.0 



r[pm] 

(1) D = 20 pm, Ca = 2.0 


Fig. 7 Radial distribution of relative platelet concentration c^ef for different tube diameters D (varying row-wise) and capillary 
numbers Ca (varying column-wise). The platelet distributions are normalised by the cross-sectional area of each bin such that 
a homogeneous distribution would lead to a constant line at Crei = 1 (denoted by horizontal dashed line). The bin size is 0.5 pm. 
Data is obtained from platelet centres and averaged between 3A^^/4 and A^. The cell-free layer thickness is shown as vertical 
dotted line. 
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• 10-3 



r[[im] 

(a) radial average velocity {D = 20 ^tm) 



_I_ I _I_ I_ 

2 4 6 8 10 

r[^m] 

(b) radial velocity fluctuation {D = 20 |um) 



(c) radial average velocity {D = 30 |am) 



(d) radial velocity fluctuation (D = 30 |um) 


Fig. 8 (a,c) Average radial platelet velocity (v^) and (b,d) its standard deviation as function of radial position r. Data is 
shown for (a,b) D = 20 pm and (c,d) D = 30 pm. Data is time-averaged after the first quarter of the simulation (i.e. between 50 
and 200 advection times), and velocities are normalised by rKBC7w The grey-shaded area denotes the CFL vicinity (2^cfl)- 
Lines are guides for the eyes. 


One can see that the gradient of at the edge of the 
CFL is larger for higher Ca. This suggests that escaping 
the CFL becomes increasingly difficult with growing Ca 
because it is less likely that a platelet is moving fast 
enough towards the tube axis before it is pushed back. 
Fig • l^can therefore explain why more platelets are seen 
leaving the CFL again in Fig. when Ca is small. 

Furthermore, the radial velocity fluctuations in 
Fig.ll in the RBC-rich region increase with Ca. This 
may result in a larger platelet diffusivity that further 
increases the lateral mobility of platelets when Ca is 
large. 

Concluding this section, one can say that non- 
diffusive platelet margination in tubes with D > 15 pm 
is only observed when Ca > 0.2. This is in good agree¬ 
ment with previous experimental results reporting that 
margination is only relevant when the wall shear rate is 
larger than 200 s“^ (Tilles and Eckstein 1987 Eckstein 


et al 1988 Bilsker et al 1989). According to the defini¬ 
tion of Ca in Eq. ([^, Ca = 0.2 corresponds to ~ 250 s“C 
This threshold value is close to the point where RBCs 
near the CEL edge are starting to tank-tread. The on¬ 
set of tank-treading therefore seems to be necessary for 


margination. The situation is different for D = 10 pm 
where platelets move to the tube wall very fast, even 
at = 0.1. A possible explanation is that many RBCs 
in Eig. |3(a)| fill the entire cross-section of the tube so 
that platelets can easily slide between them towards the 
CEL. 


4.3 Tumbling and sliding in the cell-free layer 


Once platelets are caught in the vicinity of the CEL, 
they show two main dynamical states: i) tumbling in 
the local shear flow and ii) sliding parallel to the flow 
axis (Vahidkhah et al 2014[ ). The central observable is 
the angular velocity of a platelet in a plane defined by 
the tube axis and the radial platelet position: Eig.[^ 

shows the time evolution of COr for two different platelets 
in the same simulation near the CEL. 


A freely moving platelet with axis aspect ratio p 
in a shear flow with shear rate f would be in a pure 
tumbling state leading to Jeffery orbits with average 
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An important parameter is the average relative tum¬ 
bling rate of a platelet in the CFL: 


r - ^ 

^tumb — - 

fflj 


1 I f{t)COr{t)dt 
®J 


( 12 ) 


To consider only platelets in proximity of the CFL, the 
function f{t) is defined as 



r{t) > 7?-2^cfl 
otherwise 


(13) 


Fig. 9 Angular velocity COr normalised by wall shear rate 7w 
of two platelets (grey and black) near the wall for D = 20 pm 
and Ca = 0.6. The dashed line indicates the expected average 
Jeffery frequency G)j for an unconfined platelet with aspect 
ratio 3.6. 



Fig. 10 Normalised tumbling rate Jtumb averaged over all 
platelets near the cell-free layer. Error bars indicate ensemble 
variances. Lines are guides for the eye. Note that the Ca-axis 
is logarithmic and that symbols have been slightly shifted 
along the Ca-axis to avoid overlap of error bars. 


angular velocity COj (Jeffery 1922): 




1 

= T+-. 
P 


( 11 ) 


In the present case, p = 3.6, and therefore cdj = 0.26 7 as 
indicated by the dashed line in Fig.[^ Sliding is caused 
by confinement due to the presence of the wall on one 
side of the platelet and the RBCs on the other. Only 
occasionally a platelet has “enough space” to tumble. 

As visible in Fig. tumbling events are defined 
by localised peaks in ( 0 ^, while sliding is accompanied 
by low-amplitude fluctuations of COr with zero average. 
Each tumbling event leads to a 180°- or ;r-rotation of 
the platelet, called a flip. In fact, the time integral of COr 
over one flip turns out to be close to n. A full rotation 
would therefore require two flips. 


Therefore, the integral in the denominator in Eq. (12) 
is the total time a platelet spends near the CEL: Tcfl- 
Eig.pr] shows the tumbling rate Itumb averaged over 
all platelets near the CEL. The most striking observa¬ 
tion is that tumbling is much less frequent than for 
a freely moving platelet. The tumbling rate is rather 
reduced to 10-25% of the value expected for a free 
platelet, indicating that sliding is the most likely dy¬ 
namical state in the CEL for the investigated param¬ 
eters. Vahidkhah et al (2014) reported that the slid¬ 


ing probability strongly increases with decreasing CEL 
thickness. Below ^cfl ~ 4 pm, which is the case here, 
sliding starts to occur until it eventually becomes the 
dominating dynamical mode. 

Erom visual simulation data it follows that platelet 
flips are correlated with “overtaking” events of individ¬ 
ual RBCs that move just at the edge of the CEL. On 
the one hand, a platelet is dragged along by a passing 
RBC. On the other hand, platelets And more space in 
the gaps between RBCs once an RBC has passed. This 
makes it easier for a platelet to flip, but not all RBC 
overtaking events lead to a platelet flip. 

Interestingly, Itumb is rather independent of Ca as 
seen in Eig. The tube diameter has a stronger ef¬ 
fect; larger tubes lead to a higher tumbling rate. This 
is probably related to an increase of the CEL thickness 
with D: the thicker the CEL, the smaller the confine¬ 
ment felt by the platelets. Elipping should therefore be 
easier and more frequent with increasing D but other¬ 
wise constant flow parameters (Hh? Ca). 

It is worth mentioning that the errors in Eig. 
defined by the variance of the ensemble average of all 
platelets, are relatively large. The reason is that the to¬ 
tal number of flips observed is small, and platelets show 
a wide distribution of number of flips. On average, a 
platelet flips less than ten times during one simulation. 


5 Conclusions 

Tube flow of red blood cells (RBCs) and platelets 
was simulated using a combination of the lattice-Boltz- 
mann, immersed boundary and finite element methods. 
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The aim of this work was to study platelet margination 
and subsequent platelet dynamics in the cell-free layer 
(CFL) for different tube diameters D and RBC capil¬ 
lary numbers Ca. Tubes with diameters D = 10, 15, 20 
and 30 pm were simulated, and the investigated capil¬ 
lary numbers range from C(2 = 0.1, where RBCs are rel¬ 
atively rigid, to Ca = 2.0, where RBCs are strongly de¬ 
formed and elongated. The tube haematocrit is ^37%. 

One main result is that platelet margination is facil¬ 
itated by a non-diffusive radial transport near the CFL 
edge when Ca > 0.2. This is in accordance with previ¬ 
ous experimental observations. For smaller Ca^ however, 
margination is less effective, and platelets are not irre¬ 
versibly trapped in the CFL. 

The non-diffusive platelet transport presents itself 
as a non-zero average platelet velocity in outward di¬ 
rection near the CFL edge. Platelets within the outer¬ 
most 5 pm in the tube are affected by this drift. This 
explains why platelet margination is very fast for tubes 
with D = 10 and 15 pm where essentially all platelets are 
in wall vicinity. Once a platelet has reached the CFL 
edge through diffusive motion, it is transported into 
the CFL due to collisions with tank-treading RBCs. At 
least for Ca > 0.2, this process is expected to continue 
until no platelets are left in the RBC-rich region. 

The simulations also show that the radial platelet 
velocity fluctuations are strongly suppressed in the 
CFL. The difference between these fluctuations within 
the RBC-rich region and the CFL increases with Ca^ 
which explains why platelets are less likely to escape 
the CFL with increasing Ca. 

Furthermore it is demonstrated that the predomi¬ 
nant dynamical state of platelets in the CFL is slid¬ 
ing rather than tumbling, which is in line with earlier 
findings. Due to the relatively small CFL thickness at 
~ 37% haematocrit, platelets are strongly confined, and 
tumbling events are rare. They occur 4-10 times less 
frequently as expected for an unconfined platelet. The 
tumbling rate is only a weak function of Ca^ but grows 
with D and therefore the CFL thickness. 

The presented results are in accordance with pre¬ 
vious experiments and simulations. As such, this work 
sheds more light on the platelet margination mechnism 
which is still not well understood. It is demonstrated 
how the explicit modelling of deformable RBCs is nec¬ 
essary to observe margination as an emergent effect. 
This research is thus expected to stimulate further dis¬ 
cussions and contribute to a development of a predictive 
continuum margination theory. 
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